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Abstract 

The k ■ p is a versatile technique that describes the semiconductor band structure in the vicinity of 
the bandgap. The technique can be extended to full Brillouin zone by including more coupled bands into 
consideration. For completeness, a detailed formulation is provided where the associated k-p parameters 
are extracted from the local empirical pseudopotential method in the form of band edge energies and 
generalized momentum matrix elements. We demonstrate the systematic improvement of the technique 
with the proper choice of the band edge states for the group-IV indirect bandgap semiconductors: Si, 
Ge, diamond and SiC of the 3C cubic phase. The full zone agreement is observed to span an energy 
window of more than 20 eV for Si, and 40 eV for the diamond with the 15-band pseudopotential-based 
k ■ p approach. 
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1. Introduction 

For the electronic and optical processes in semiconductors involving states near the bandgap, the k ■ p 
technique has been the first resort for many researchers [1]. The k ■ p method was introduced by Bardeen [2] 
and Seitz [3] . Luttinger and Kohn extended the technique to degenerate bands to govern the valence band 
edge structure of common group-IV elemental and III-V compound semiconductors as described by a 3 x 3 
matrix [4]. In the case of direct bandgap semiconductors, appending the lowest conduction band as well as 
including the spin-orbit split-off within the valence band lead to a 4 x 4 matrix which is also known as the 
Kane's model [5]. Soon afterwards Pidgeon and Brown extended this to 8 x 8 matrix by including the spin 
for all bands in the presence of an external magnetic field [6] . These measures have enlarged the number of 
coupled bands and also improved the overall agreement of the band effective masses with the experimental 
values. On the other hand, their validity spanned only about 15% of the full Brillouin zone [7] which can be 
for instance sufficient for the description of not too small excitons. However, they fall short for applications 
involving the larger portion of the Brillouin zone such as the optical absorption spectra or the high-field 
transport characteristics. The classical work in the direction to extend the multiband k ■ p to full zone is 
due to Cardona and Pollak [8] which resulted in a 15 x 15 k ■ p Hamiltonian. After this period of progress 
over the technique, researchers especially for device physics applications have routinely employed the 4-band 
and the 8-band variants of k ■ p together with the envelope function approximation for incorporating the 
confinement effects [9]. Wood and Zunger have termed this combination as the "standard model" and they 
have compared it with the empirical pseudopotential method (EPM) [10]. Both qualitative and significant 
quantitative errors were identified mainly due to L— and AT— derived states which were not accounted by 
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the 8-band k ■ p approach [10]. The particular concerns on the use of the envelope functions were relieved 
by Foreman [11] as well as by Burt who also clarified the correct operator orderings [12]. 

For the nanotechnology applications involving quantum dots and nanocrystals, the device sizes are re- 
duced in all three dimensions which results in the sampling of the full Brillouin zone. Once again there is 
the necessity to extend the multiband k ■ p to full zone for technologically-important semiconductors. The 
increase of bands also results in the proliferation of the so-called Luttinger-like parameters which cannot be 
easily extracted from experiments as was the fortunate case for 3-bands where effective masses were directly 
related [1]. One possibility is to make use of the tight-binding model as done by Jancu et al. who have 
systematically determined the 14-band k ■ p parameters from the 40-band tight-binding band structure [13]. 
Another route is to resort to the classical Cardona-Pollak methodology [8] where the empty lattice plane 
wave basis is used; Richard et al. followed this approach which lead to 30-bands taking into account the 
spin-orbit coupling. An even simpler alternative is to make use of the already existing EPM band structure 
data to generate the multiband k ■ p parameters [14]. It is this last choice that we pursue in this work to 
obtain the multiband k ■ p band structure for group IV indirect bandgap semiconductors: Si, Ge, diamond 
and SiC of the 3C cubic form. Our basic motivation behind this choice is that these semiconductors act 
as embedded core materials for nanocrystals where the confinement effects alleviate the poor luminescence 
properties of their bulk counterparts [15]. However, the bulk electronic structure is nevertheless inhereted to 
low-dimensions, with the further complication from the interface states. We believe that for certain nanotech- 
nology applications, full zone k ■ p band structure that can easily be extracted from EPM as demonstrated 
in the following sections may still be a suitable choice compared to more demanding atomistic approaches. 



2. Theory 



The bulk crystal Hamiltonian described by a local one-electron effective crystal potential, V^- ta j(r) is 
given by 



2toq 



+ WO 



(1) 



where p = —ihV, mo is the free-electron mass, n is the band index, and k is the crystal wavevector. 
Inserting for the Bloch wave functions ^P n j:(r) — ~^ elkr u nf{f)i where V is the crystal volume, we obtain 
the well-known the k ■ p expression 
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where H = \^ + V^ ta i(r). Note that the cell-periodic functions {u n ^(r}} satisfy a Sturm-Liouville type 
eigenvalue equation with the eigenfunctions being complete when all (i.e., infinite) bands are included at 
a fixed wavevector k, say k - Based on this completeness at a chosen fc we can expand an arbitrary 
cell-periodic function, u.^{r) as 



Aft 



r(r) = b r r u r (r) ■ 



(3) 



m— 1 



Inserting this form into the bulk Hamiltonian given by Eq. (1) leads to the following eigenvalue equation 
that becomes exact as Nt, — > oo: 



where 
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corresponds to the momentum matrix element. 

For computations, obviously one cannot include infinite bands, so it is more viable to build up an 
expansion set from a number of band and wavevector combinations chosen with a physical insight mainly 
reflecting the band extrema states. In other words, rather than using a fixed ko we form an expansion over 
various band and wavevector indices, {to, ki} as 



{m,ki} 



(7) 



Since the cell-periodic functions are not necessarily orthonormal at different wave vectors, we define their 
overlap as 

(8) 
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Hence, the expansion coefficients & J -jf m jf j satisfy the following equation 
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Projecting to \u n j; ) results in a generalized eigenvalue equation for the expansion coefficients 
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(10) 



where {n, k s } also belong to the same chosen expansion set as {to, ki}; here 
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with 
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(12) 



corresponding to the generalized momentum matrix element linking indirect transitions. 

In summary, one solves Eq. (10) for the desired band j and the wavevector k with the only ingredients 
being the band energies, E m (ki) and the generalized momentum matrix elements, p n j: m £. . These can easily 
be extracted from a local pseudopotential approach; the cell-periodic functions can be expressed in the 
Fourier representation as 

-„*(-) = 4=E^(G)e^. (13) 



G 

where {G} correspond to reciprocal lattice vectors and ilo is the volume of the primitive cell. Hence the 
overlap (Eq. (8)) and the generalized momentum matrix elements (Eq. (12)) arc given in terms of Fourier 
coefficients -B m g. (G) of the cell-periodic functions as 



A r t = VB*. (G)B t (G) , 
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and 
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respectively. 
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Figure 1. (Originals in color) EPM (black lines) vs. 4-band k -p (red symbols) band structure for Si. 



3. Results 



We apply this basic idea to group IV indirect bandgap semiconductors: Si, Ge, diamond and SiC of the 
3C cubic phase. Their local EPM form factors are taken from Refs. [16, 17] for Si and Ge, respectively and 
from Ref. [18] for the diamond and SiC. Spin-orbit coupling effects are neglected which can actually become 
significant for the case of Ge [17]. Another important technical remark is about the EPM cut off energies: 
we observed that even though the EPM band energies (i.e., eigenvalues) converge reasonably well with cut 
off energies as low as 5-10 Ry, the corresponding Bloch functions (i.e., eigenvectors) require substantially 
higher values to converge. The results to follow are obtained using 14, 16, 33, and 22 Ry for Si, Ge, diamond 
and SiC, respectively. 

Using silicon as a testbed, we demonstrate how to choose and improve the k ■ p basis set. First, we start 
with the four-band k ■ p which is known to yield reasonable results for direct bandgap semiconductors like 
GaAs. These four states are taken to be the highest three valence band (VB) states at the F point and the 
lowest conduction band (CB) state at one of the X points, (1,0,0). The resultant band structure displayed 
in Figure 1 indicates that apart from the four band edge energies even the band curvatures turn out to be 
wrong. To improve this situation, the eight-band k ■ p is generated by employing the highest four VB and 
the lowest four CB states, all from the V point. The agreement with the EPM band structure shown in 
Figure 2 is once again not acceptable due to remarkable deviation away from the T point. Also note that 
the heavy hole band acquires the wrong curvature along both the T — L and T — K directions. To remedy 
these shortcomings, we form a set by including from the T point the band indices 1 to 8 (as in the previous 
case), from the X point bands 3 to 6, and from the L point bands 3 and 4, and finally from the K point the 
fifth band making altogether 15 states. Here the band index 4 corresponds to highest VB, and 5 to lowest 
CB. The band structure of this pseudopotential-based 15-band k ■ p approach is shown in Figure 3 where 
excellent agreement is observed over an energy window of about 20 eV. 

Using the same states for the 15-band k ■ p approach for the other indirect bandgap semiconductors Ge, 
diamond and SiC, we reach the same performance with Si as displayed in Figs. 4-6, respectively. For the 
diamond which is a wide bandgap semiconductor, the full zone agreement is observed to span an energy 
window of more than 40 eV. 
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Figure 4. (Originals in color) EPM (black lines) vs. 15-band k ■ p (red symbols) band structure for Ge. 




Figure 5. (Originals in color) EPM (black lines) vs. 15-band k ■ p (red symbols) band structure for diamond. 
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Figure 6. (Originals in color) EPM (black lines) vs. 15-band k ■ p (red symbols) band structure for SiC in 3C cubic 
phase. 

4. Conclusions 

The multiband k ■ p technique can be systematically improved by including more bands. This is particu- 
larly needed for indirect bandgap semiconductors where 4- and 8-band approaches fail. The band coupling 
parameters which are in the form of generalized momentum matrix elements and the associated band edge 
energies for the selected states can be extracted from EPM. The technique can also benefit from ab initio lo- 
cal pseudopotcntial band structure provided that the bandgap shrinkage due to local density approximation 
is corrected [19]. We believe that for certain nanotechnology applications, full zone k ■ p band structure that 
can easily be generated with the inputs from EPM as demonstrated in this work may become the suitable 
choice compared to more demanding atomistic approaches. Therefore, what remains to be done is to check 
the performance of the proposed 15-band k ■ p framework on several low-dimensional applications and to see 
whether the "farsightedness" of the conventional k-p approach with the small number of bands [20] is cured. 
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